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Femtoscopic measurements at RHIC have been hailed as a source of 
insight into the bulk properties of QCD matter. However, hydrodynamic 
models, which have been successful in reproducing other observables have 
failed to satisfactorily explain femtoscopic two-particle correlation mea- 
surements. This failure has been labeled the "HBT puzzle". In this talk, 
I present correlations using a hybrid model composed of a viscous hydro- 
dynamic module and a hadronic cascade. After incorporating early accel- 
eration, viscosity and a stiffer equation of state, the extracted source sizes 
come much closer to experiment. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.70.Pq 

At the dawn of the RHIC era, the most sophisticated, and seemingly 
most realistic, transport models failed miserably to match source sizes as 
inferred from experimental two-particle correlations measurements. These 
models [H [2j [3] applied a hydrodynamic code to model the super-hadronic 
stage, which was then coupled to a hadronic cascade whose purpose was 
to simulate the low density stage and breakup. Seemingly less realistic 
cascade models were not completely satisfactory, but performed significantly 
better [U [5J El \7\. These models ignored the partonic phase altogether, 
or at least ignored the softness of the equation of state associated with 
the phase transition. The failure of the hydrodynamic models picked up 
the name "HBT puzzle" (HBT refers to Hanbury-Brown and Twiss who 
developed interferometric source size measurements with photons [8j). A 
second feature of the puzzle concerned parametric fits, such as the blast- 
wave model, which could come close to the data but only with unphysically 
high breakup densities. Inspired by this puzzle, we have developed a viscous 
hydrodynamic model coupled to a hadronic cascade [9], and have studied 
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the effect of several improvements or modifications of the initial treatments, 
and found that the femtoscopic data from RHIC can be well reproduced 
without either invoking any contentious ideas, or setting parameters outside 
usual expected values. Considered individually, none of these modifications 
explained more than half the original discrepancy. Instead several changes 
conspired to push source sizes in the same direction, mainly by making the 
reaction more explosive which leads to smaller source sizes. 

Since many talks from this workshop have already described how fem- 
toscopic source sizes can be extracted from data, and how they are theo- 
retically related to the emission history of a heavy ion collision, I shall be 
brief, and mention a few key points. First, the three dimensions R ou t,R s ide 
and -Riong refer to the outward (parallel to the momentum of the particle), 
sideward (perpendicular to the beam and to the particle's momentum) and 
longitudinal (along the beam axis) directions. The sizes represent Gaus- 
sian fits to the spatial size and shape of the outgoing phase space cloud of 
particle's of a specific momentum. Since the initial fireball has little ex- 
tent along the beam direction, the size i?i on g can be related to the time at 
which particles are created. If the emission is sudden, the two transverse 
dimensions tend to be similar, but for a longer duration of emission, par- 
ticles emitted earlier get ahead and i? out becomes larger than i? s ide- For 
central collisions, and for measurements at mid-rapidity, each dimension is 
a function of pt only. The aforementioned hydrodynamic models tended to 
over-predict i?i ong and the -R ut/-R s ide ratio, suggesting that the explosion 
was more rapid, and with a more sudden emission than that portrayed by 
the models. 

The principle changes to the model responsible for improving the com- 
parison with data are: pre-equilibrium flow, using a stiffer equation of state, 
and incorporating viscosity. A modest improvement also entails if one more 
realistically models the tt—tt interaction. Additionally, we affirm the findings 
of [12] and see that the shape of the initial profile can significantly affect 
the source sizes. In order to illustrate the effects, Fig. [1] first presents a 
benchmark calculation, which has no pre-equilibrium flow, no viscosity and 
too stiff of an equation of state. The source sizes poorly match the data, and 
are similar to other models with the same lack of features. The -R ut/-Rsidc 
ratio is approximately 40% too high, and i?i ong is also over-predicted. By 
moving the start time of the calculation from 1.0 fm/c to 0.1 fm/c, the colli- 
sion becomes more explosive, which shortens the extent and duration of the 
pion emission, and leads to source sizes significantly closer to the data as 
shown in Fig. [TJ Similar effects have been seen in similar contexts in both 
hydrodynamic and microscopic models |X3|. El [12] . Pre-equilibrium flow is 
inevitable, as flow derives from the conservation of the stress-energy tensor, 
and its strength is largely independent of equilibration or the microscopic 
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Fig. 1. Gaussian radii reflecting spatial sizes of outgoing phase space distributions 
in three directions: R on t, Rside and i?i on g- Data from the STAR collaboration (red 
stars) are poorly fit by a model with a first-order phase transition, no pre-thermal 
flow, and no viscosity (solid black squares). Correcting for all those deficiencies, and 
using a more appropriate treatment of the relative wave function in the Koonin 
equation [THl [H] brings calculations close to the data (filled black circles). The 
sequential effects of including pre-thermal acceleration (open blue squares), using 
a more realistic equation of state (open green diamonds), and adding viscosity 
(open cyan triangles) all make substantial improvements to fitting the data. An 
improved relative wave function yielded modest improvements (compare open cyan 
triangles to filled black circles). 

state of matter during the first fm/c |14j . 

The second change we consider is improving the equation of state. The 
equation of state used in the benchmark calculation was first order, with 
a large latent heat, L = 1.6 GeV/fm 3 , similar in spirit to what was used 
in previous studies [H [U [3]. During the mixed phase, the speed of sound 
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is zero and the pressure and temperature stay constant, which would lead 
to extended lifetimes, especially if the initially energy density were lower, 
and close to the maximum for the mixed phase |X5|, [16] . Lattice calculations 
suggest a noticeably stiffer equation of state with no first-order transition. 
Figure Q] shows significantly improved source sizes from a calculation where 
the soft region is cut in half and the speed of sound, c 2 = dP/de, is set to 
0.1c. As expected, the stiffer equation of state increases the explosivity and 
results in improved source sizes as illustrated in Fig. [TJ 

The third modification to the hydrodynamic model is to include vis- 
cosity. The Navier-Stokes equation gives the following form for the stress- 
energy tensor in the fluid rest frame, 

Tij = P5 i:j - r, {d lVj + djVi - (2/3)V • v<%} - CV • v, (1) 

where r\ and £ are the shear and bulk viscosity coefficients respectively. In 
the early stages of a RHIC collision the velocity gradient is far stronger along 
the z axis, and the shear correction is negative for T zz and positive for the 
transverse components of pressure, T xx and T yy (The correction for shear is 
traceless). This increases transverse acceleration which again significantly 
improves agreement with experimental source sizes as shown in Fig. [H and 
was also seen in |17] , The shear viscosity in the partonic phase used here 
is twice the KSS limit [18]. A bulk viscosity was also applied that peaks 
in the mixed phase, qualitatively consistent with both phenomenological 
models and with lattice extrapolations [19j [20], but didn't lead to strong 
modifications to the source sizes. More details describing how bulk and 
shear viscosity are implemented can be found in [9]. 

Finally, we apply an improved treatment of the ir — ir interaction used in 
calculating the relative wave function in the Koonin formula |10[ [TT] used 
to generate correlation functions, 

C(k*,q) = J d 3 r 5(k t ,r)|0(q,r| 2 (2) 
S(k t ,r) = 

Um J d^r 2 f(k' t , r[,t')f(k' t , - n' + r 2 Q 

t'^oo Jd 3 nd 3 r 2 /(k t ,ri,t)/(k t ,r 2 ,t) 

Here, f(k t ,r, t) represents the phase space density for pions of the rele- 
vant momentum, and the source function S(bfk t , r) provides the normalized 
probability for two particles of the same momentum, kt, to be asymptoti- 
cally separated by r in the pair rest frame (The primes refer to quantities in 
that frame). For the calculations presented thus far, a simple symmetrized 
plane wave for was used for the relative wavefunction, \(p\ 2 = 1 + cos(q ■ r), 
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which makes calculation of the correlation functions simple. More impor- 
tantly, the correlation functions for Gaussian sources are then Gaussians in 
relative momentum, which makes the search for the best fit easy. If the 
sources were truly Gaussian, the radii extracted from this method would 
match those from a procedure using the full relative wave function including 
the strong and Coulomb interaction between pions. For the final calculation 
portrayed in Fig. Q]more realistic correlation functions were generated using 
the full wave function. These correlations were then fit to Gaussian sources 
using the Bowler-Sinyukov method |2H [22] used in experimental analyses 
to mimic Coulomb effects. This results in a slightly improved fit to the 
experimental source sizes in Fig. [TJ 

After the four improvements listed above, the model now reproduces the 
experimental source sizes to better than 10% for all radii, which would seem 
to be within the systematic error quoted by the experiments, and close to 
the systematic uncertainty associated with the basic phenomenology. Thus, 
no single effect accounted for more than half of the HBT puzzle. Instead, 
the explanation amounted to a conspiracy of several effects, each of which 
increase the explosivity, and pushed the source sizes toward the data. An- 
other modification, which will not be illustrated here, is to change the shape 
of the initial energy density profile. More compact profiles make the colli- 
sion more explosive |12j . The profile applied here is that of the wounded 
nucleon model [23], but color-glass profiles 

The second half of the HBT puzzle concerned the breakup densities in- 
ferred from the parameterized models \14\ [25] . Over a thousand hadrons are 
produced in a central collision. For a blast-wave femtoscopic fitting exer- 
cises typically suggested radii near 12 fm, and a breakup time near 9 fm/c. 
The breakup volume would then be ttR 2 t, or ~ 4000 fm 3 , corresponding 
to densities near 0.25 hadrons per fm 3 , almost twice nuclear density. Given 
that hadronic cross sections tend to be ~ 30 mb, or 3 fm 2 , the mean free 
path would be below 2 fm, which is far below the overall system size. The 
solution to this puzzle can be understood by viewing Fig. [2j which displays 
points sampling the final emission points from they hyrodynamic/cascade 
model for pions whose asymptotic momenta is 300 MeV/c pointed along the 
x axis. The x and t coordinates are clearly positively correlated, as expected 
for emission from an expanding surface, and are similar to what was seen 
in AMPT [7|. Since the source sizes characterize the shape of the outgoing 
phase space distribution for particles of this particular velocity, this correla- 
tion allows later produce pions to be produced in proximity to their earlier 
produced counterparts. If the particle emission had been confined to the 
neighborhood of a line with the slope of the velocity (see dashed lines), the 
emission duration would have been irrelevant, and R ut/Rside would have 
been even lower than found here. The positive x — t correlation allows the 
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Fig. 2. Final emission positions and times for particles with transverse momentum 
of 300 MeV/c along the x axis. Emission has a positive correlation between position 
and time, though lags behind the slope of the velocity (illustrated by dashed lines). 
Due to the positive x — t correlation, emissions of longer duration can still result in 
phase space clouds that are compact along the outward direction, with i? ut / -Rsidc ~ 
1. 



-Rout / Rs\de ratio to be near unity, even though the average emission time was 
near 20 fm/c. The parameterized forms do not properly account for this 
correlation, hence they provide a misleading picture of the breakup process. 

The success of the fits in Fig. Q] are encouraging, but they immediately 
raise the question of whether the changes to the model destroyed the ability 
to fit spectra. To that end, Table [1] compares the mean pt values for the 
final model to experimental values from STAR and PHENIX. For protons, 
pions and kaons, the models value lied between the values reported by the 
two collaborations. Given that the previous hydrodynamic treatments also 
fit these values, it is surprising that the modifications presented here, all of 
which made the collision more explosive, did not result in more explosive 
spectra. In particular, we had expected the mean pt for protons to be over- 
predicted. Some of this may be due to the addition of bulk viscosity, which 
had little affect on HBT, but did lower the mean pt. It is also possible, this 
agreement will disappear once the spectra are compared more carefully. In 
particular, it should be mentioned that the mean pt values from the model 
included products of weak decays. Filtering such products in a manner 
consistent with experiment might destroy the agreement. 

Along with spectra and femtoscopic correlations, the third class of hadronic 
bulk observables are large-scale correlations, which are associated with col- 
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p,n,p,n 


STAR 


26 




422 ± 22 


719 ± 74 


1100 ± 110 


PHENIX 


27 




453 ± 33 


674 ± 78 


954 ± 85 


Hydro+Cascade 


433 


714 


1027 



Table 1. The mean (p t ) in MeV/c for central collisions for pions, kaons and pro- 
tons. Only charged species were used in the PHENIX analysis, and only negative 
hadrons were used for STAR. Since the original energy density was chosen to fit 
the multiplicity, the success seen here in reproducing the mean p t suggests that the 
spectral shapes are probably well described. 

lective flow. The elliptic-flow variable should be enhanced by the inclu- 
sion of pre-equilibrium flow, and by using a stiffer equation of state. Shear 
viscosity would be expected to lower the V2, so it would not be surpris- 
ing if these modifications resulted in little net change. Unfortunately, the 
hydrodynamic code used here assume radial symmetry, which precludes el- 
liptic flow analysis. Nonetheless, the result found here provides significant 
hope that the entirety of soft hadronic observables can be reproduced with 
a single model. Doing so would represent a significant milestone for the 
field. This would be no means validate such a model, or its parameters, but 
it would demonstrate that a rigorous statistical assault on the data would 
result in a non-null region. 
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